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Aims. We present the first results of a wide field survey for cool brown dwarfs with the MegaCam camera on the CFHT telescope, 
the Canada-France Brown Dwarf Survey, hereafter CFBDS. Our objectives are to find ultracool brown dwarfs and to constrain the 
■ field-brown dwarf mass function thanks to a larger sample of L and T dwarfs. 

Methods. We identify candidates in CFHT/MegaCam i' and z' images using optimised psf-fitting within Source Extractor, and follow 
them up with pointed near-infrared imaging on several telescopes. 

Results. We have so far analysed over 350 square degrees and found 770 brown dwarf candidates brighter than z^ fl =22.5. We currently 
have ./-band photometry for 220 of these candidates, which confirms 37% as potential L or T dwarfs. Some are among the reddest and 
— ' farthest brown dwarfs currently known, including an independent identification of the recently published ULAS J003402.77-005206.7 

| and the discovery of a second brown dwarf later than T8, CFBDS J0059 10. 83-01 1401.3 . Infrared spectra of three T dwarf candidates 

, confirm their nature, and validate the selection process. 

Conclusions. The completed survey will discover ~100 T dwarfs and ~500 L dwarfs or M dwarfs later than M8, approximately 
doubling the number of currently known brown dwarfs. The resulting sample will have a very well-defined selection function, and 
will therefore produce a very clean luminosity function 
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120071) . Follow up of these discoveries, and of the fewer brow n 
dwarfs identified as companions to stars (e.g. lScholz et al.l2 003). 
has lead to spectacular advances on (1) interaction between mat- 
ter and radiation in cool dense and complex turbulent atmo- 
sphere, where molecules and dust form and dissipate; (2) stellar 
and planetary formation and (3) galactic structure, wit h the first 
chara cterization of the substellar mass function (e.g. [Chabrier 
d200ll) . for an extensive review). 

Much interesting work however remains to be done, and 
the advent of wide field cameras on large telescopes makes 
an unprecedented volume of the Milky Way accessible for 
brown dwarf searches. Here we use two large surveys with 
MegaCam 1 on the CFHT telescope, the Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS 2 ) and the Red-sequence 
Cluster Survey 2 (RCS-2; lYee et al.1 (|2003)), and complement 
them by additional observations to address three areas of brown 
dwarf physics: 

(1) Detection and physics of ultracool brown dwarfs (T eff < 
1000 K). As of today, observed stellar and substellar atmo- 
spheres cover a continuum of physical conditions from the 
hottest stars (~ 100 000 K) to the coolest known brown dwarf 
dDelorme et al.l 120081 ~ 625 K,). There remains a sizeable 
temperature gap, between these coolest brown dwarfs and the 
~100 K giant planets of the Solar System. Besides their intrinsic 
interest, ultracool brown dwarfs provide analogs to these planets 
that are not encumbered by the glare of a bright star. This will 
greatly help guiding the design of direct planet detection exper- 
iment, which currently have to rely on unvalidated models. 

(2) Brown dwarfs in the thick disc and spheroid. The LSR 
(Lepine-Shara-Rich) and 2 MASS surveys have recently iden- 
tified the first L subdwarfs (Burgasseret al. 2003; Lepin eet al.l 
120031) . i.e. low metallicity brown dwarfs from the galactic halo 
population, from samples of a few hundred L dwarfs. The 
MegaCam survey will reach further down the halo luminosity 
function and may find a few T-type subdwarfs. 

(3) Statistics of brown dwarfs of intermediate temper- 
ature (1000 to 1500 K). Current estimates of the substellar 
Galactic mass function suggest that in the disk of the Galaxy 
brown dwarfs a re about as numerous as stars (e.g Chabrier 2001; 
ICruz et al I l2007h . That mass function however has significant 
statistical noise, which reduces its power as a constraint on star 
formation and galactic structure theories. At these low effective 
temperatures the final luminosity functions from the DENIS, 
2MASS and SDSS surveys will all retain significant Poisson 
noise The samples from individual searches are not easily com- 
bined since they are affected by different selection biases, so 
only a fraction can b e used to define a robust luminosity function 
(e.g. lCruz et al.l2007h . By almost doubling the number of known 
brown dwarfs, from a single survey with a well understood se- 
lection function, we will provide significantly tighter constraints 
on the luminosity function. 

The present paper describes the overall strategy of our 
brown dwarf search. Sec 2 discusses the observational proper- 
ties which we use to identify these extremely red objects, and 
presents the corresponding observational material, while Sec. 3 
describes how we generate a candidate list with minimal con- 
tamination from both observational and astrophysical artefacts. 
Sec. 4 describes the characteristics of the resulting candidates 
and presents spectra for a few T dwarfs identified early-on. We 
conclude with a discussion of the expected results for the com- 
pleted survey. 



1 http://www.cfht.hawaii.edu/Instruments/Imaging/MegaPrime/ 

2 http://www.cfht.hawaii.edu/Science/CFHTLS/ 



2. Observations 

2.1. Observational properties of brown dwarfs 

Field brown dwarfs are extremely cool objects, with a temper- 
ature range whic h currently extends from ~ 2500K (early L) to 
~ 625K (late D dDelorme et al.ll2008tlGolimowski et ai1l2004t 
IWarren et al.l l2.Q07b . Even cooler, yet to found, brown dwarfs 
should close the temperature gap between late type T dwarfs 
and solar system giant planets(~ 100K). Brown dwarfs spectra 
very much differ from a black body, and have considerable 
structure from deep absorption lines and bands. Their spectral 
energy distribution (in vF v units) peaks in the near infrared 
(hereafter NIR), particularly in the J photometric band, and they 
are most easily detected in that wavelength range. Their pure 
NIR JHK colours however do not very effectively distinguish 
them from other classes at modest S/N ratio. Brown dwarfs are 
more easily recognised by including at least one photometric 
band below 1 /mi, since their steep spectral slope at those 
wavelength produce very distinctively red i' - z' and z' - J 
colours. As one re cent example, the T 8.5 ULAS 0034 has 
(/' - z)' AB >3.0 (5 o-. lDelorme et al.ll2008l) . and at any S/N ratio 
where it is safely detected it cannot be confused with anything, 
except a slightly earlier T dwarf or a z=6 quasar. The UKIDSS 
discovery observation however was less than 3<x away from the 
K dwarf locus. Since K dwarfs outnumber T dwarfs by orders 
of magnitude in any flux limited sample, that distance would 
have been woefully insufficient for a secure identification. The 
non-detection of ULAS 0034 at i' and z' in th e deep SPSS 
stripe 82 played a major role in its identification bv lWarren et al.l 
d2007l) ). and other near-IR searches for brown dwarfs similarly 
use some <l//m imaging to weed out their contamination. 

We t ake advantage of the wide field of view of the MegaCam 
camera (Boulade et al. 2003) on the CFHT telescope, and of the 
trove of observational material obtained with that instrument, to 
select brown dwarfs on their i' - z' colour. The i' - z' colour has 
excellent dynamics for brown dwarfs, varying from 1 .7 to 4.0 be- 
tween mid-L and late-T (Fig.[3}. It therefore provides (at least at 
high S/N ratio) a good spectral type estimator. At the high galac- 
tic latitude of our survey, the i' - z' colour distinguishes brown 
dwarfs from almost every astronomical source type, but it leaves 
one contaminant, quasars at z > 5.8. Those are of considerable 
interest in their own right, but ne ed to be distingu ished from the 
brown dwarfs. As first shown bv lFan et al.l d2001l) . the i' - z' vs 
z' — J colour/colour dia gram very effective ly separates the two 
populations (Fig. [T]and Will ott et al.l d2005)). The very red i' - z' 
of high redshift quasars is caused by deep Lyman a absorption 
on a relatively flat intrinsic spectrum), and they therefore have a 
more neutral z' - J. The spectral distribution of brown dwarfs, 
in contrast, keeps rising steeply into the J band. We therefore 
complement our MegaCam i' and z' photometry by pointed J- 
band imaging of the candidates selected on i' - z' ■ Besides pin- 
pointing the (few) quasars, the /-band photometry very effec- 
tively rejects any remaining observational artefact, as well as the 
(more numerous) moderately red stars scattered into the brown 
dwarf/quasar box by large noise excursions. 

2.2. Synthetic colours 

Each square-degree MegaCam image contains a few hundred 
thousand objects, of which at most a few are brown dwarfs. 
We thus need to strike a careful balance between sample 
completeness and contamination. To tune this compromise we 
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Fig. 1. z' - J and i' - z' colours of brown dwarfs and quasars. 
These synthetic colours were computed for the MegaCam i' and 
z' photometric system and the NTT SOFI J s i wr t system. The 
bold lines represents the average colours of dwarfs (dashed) and 
quasars (solid). Symbols marks the individual colours of known 
brown dwarfs. The thin grey lines mark the colour evolution of 
indivi dual synthetic quasars (described in detail in Willot t~e~t alj 
(2005)) as their redshift ranges from 5.5 to 6.7. The quasar red- 
shift increases for redder colours. The vertical line at i' - z'=1.7 
marks our i' - z' selection criterion. The quasar selection box is 
also marked (to the lower-right of the solid line). 



need a precise knowledge of the colours of brown dwarfs and 
quasars for the e xact instruments and filters used in our survey. 
As discussed in Wil lott et a D (120051) . these colours are known 
for some photometric systems, in particular SDSS and 2MASS, 
but the filters and quantum efficiency curves of MegaCam 
are notably different (Fig|2]i. This is particularly significant for 
brown dwarfs and quasars. Due to their highly structured spectra 
a modest change to a response curve can produce significantly 
different colours when it includes or excludes a major absorp- 
tio n band or emission line. We update the synthetic colours 
of I Willott etall (I20051) . using additional brown dwarfs spectra 
which have become available since 2005, and adding the many 
near-IR instruments and filters which we use for the /-band 
imaging. W e use publica lly availa bl e spectra (from S. L egged:' s 
website 3 , iMartinetail dl999l): iKirkpatrick et alj d2 000); 
iGeballe et all d2001l); iLeggett et all d2002l); IBurgasser et alj 
d2003l): iKnapp et alj d2004l) : iGolimowski et alj d2004l) : 
IChiuetalJ (12006)) of o ver 60 brown d warfs with spectral 



types LI to T8 (on the IBurgasser et alj d2006l) spectral typ e 
scale) and the synthetic quasar spectra of Will ott et all d2005). 
We compute th eir synthetic MegaC am i' and z' photometry in 
the AB system (Fukugit a et alj!996l) using detector quantum ef- 
ficiency and transmission curves for the atmosphere, telescope, 
camera optics, and filters (cf figur^H, obtained from the CFHT 
web page. Figure. [3]displays the resulting colours as a function 
of the spectral type. 

We similarly synthesize /-band photometry for each of the 
instruments and / filters used in the /-band follow up. These 
instruments have significantly different response curves, which 



must be taken into account to obtain homogeneous selection cri- 
teria. We found, in particular, that brown dwarfs colours which 
include / photometry obtained at the NTT with SOFI and its 
(default) wide / filter are not as red as we initially expected: 
contrary to most / filters, its wide bandpass includes water va- 
por bands which are strongly absorbed in L and (particularly) T 
dwarfs. As a result, the i' - J and z' — J colours which use this 
filter are bluer by ~ 0.15 mag for early L and ~ 0.5 mag for 
late T. After we realised this we switched our SOFI observations 
to the alternate J short filter, which better separates T dwarfs from 
quasars. We use the synthetic colours to shift our selection boxes 
according to the actual filter. 



Compared transmission in i' and z* band for CFHT and SDSS 
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Fig. 2. Compared spectral response functions of the CFHT (thick 
lines) and SDSS (thin dashed lines) instruments for their i' (dark 
blue) and z' (red) filters. These factor in the average atmospheric 
transmission of the two observatory sites, the telescope reflectiv- 
ities, the transmissions of the camera optics and filters, and the 
quantum efficiencies of the CCDs. Contrary to the SDSS band- 
passes, the CFHT i' and z' filters overlap significantly, leading to 
less contrasted colours. 



2.3. Optical data 

Our survey for brown dwarfs, the Canada-France Brown Dwarf 
Survey (CFBDS), builds upon two major MegaCam surveys, 
the Canada-France-Hawaii Telescope Legacy Survey (hereafter 
CFHTL S) and the Red- sequence Cluster Survey 2 (hereafter 
RCS-2, lYee et alJ 120071) and, where necessary, complements 
their filter coverage with the additional MegaCam observations 
needed to obtain pairs of i' and z' images. Figure H] summarizes 
the currrent sky coverage of our optical survey data. 

The CFHTLS survey has three components, named Deep, 
Wide, and Very Wide, ( described in detail on the CFHTLS web 
page, http://www.cfht.hawaii.edu/Science/CFHLS/), and we use 
all three. 

CFHTLS Deep This deepest component of the CFHTLS 
covers four high galactic latitude and low extinction 1 
square degree fields in the u'g'r'i'z' filters, and it is pri- 
marily motivated by the SNLS type la supernovae search 
dPain & SNLS Collaboration! 120031) . The expected total expo- 
sures times per field of the completed survey will be 66h at z' 



and 132h at i' with 10 cr depths of z' A 



AB 



-26.25 and i' AB 



-27.65. 
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Table 1. Characteristics of the optical surveys used by CFBDS. 
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The detection limits correspond to lOcr, as needed for 10% precision photometry. 




Fig. 4. Sky chart of sky area covered by CFBDS so far. Black curve marks the galactic plane while dotted curves mark +30° and 
-30° galactic latitude. 



This sensitivity is sufficient to identify mid-L dwarfs out to 
1300 parsecs. At such distances, and at high galactic latitude, 
the thick disk becomes dominant and we therefore have good 
hopes to probe its brown dwarf population. The present analysis 
uses the T0003 release (T0004 was very recently released, 
but we have yet to analyse and follow-up the corresponding 
detections) which has 10 cr depths of z' AB ~24.7 and i' AB ~26.25. 
Thanks to their large number of contributing exposures and their 
extended time base, the T0003 images are extremely clean, with 
essentially no contamination by cosmic rays or bad pixels, or by 
variable or moving objects such as supernovae and asteroids. 

CFHTLS Wide This component of the CFHTLS will 
cover 186 deg 2 divided between four high galactic latitude low 
extinction fields, in u*g'r'i'z' , and it is primarily motivated by 
cosmological weak lensing. We have analysed the 20 deg 2 that 
have both i' and z' coverage in the T0003 CFHTLS release (/' 
coverage is considerably more extensive, due to priorities set by 
the main drivers of the Wide survey). The average 10 cr depths 
are z' AB ~23.8 (for total exposure times of 7200s) and i' AB ~24.75 



(for total exposure times of 4300s), with small field to field 
variations due to seeing and sky background fluctuations. Of 
the three components the Wide probes the largest volume. 
The Wide images have enough coadded subexposures (9 for z' 
and 7 for i') to reject all cosmic rays and bad pixels, and the 
overall exposure times are sufficiently long to eliminate all but 
the slowest moving objects. The i' and z' images on the other 
hand are usually not contemporaneous, and variable sources (in 
practice mostly supernovae) which are serendipitously bright in 
the z' can erroneously pass our i' - z' colour filter. Those need 
to be eliminated at a later stage. 

CFHTLS Very Wide (VW) This shallowest component 
of the CFHTLS is motivated by transneptunian objects and 
was initially set to cover 1000 deg 2 in the ecliptic plane with 
g'r'i' images. It was later downsized to ~250deg 2 when it was 
realised that the three components could not all be completed 
within the allocated time. We use the ~ 150 deg 2 from the 
Very Wide with absolute value of the galactic latitude above 
30 degrees to ensure low absorption (see table Q]). The average 
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Synthetized colors of known L and T dwarfs 
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Spectral Type(L,T) 

Fig. 3. i' - z' synthetic colour versus spectral type for the 
MegaCam photometric system. The colour for a sp ectral type 
is the average over the b r own d warf spectra from Chiu et all 
(120061) . iGohmowski et alJ d2004 and iKnapp et"all (12004 with 
that spectral type, and the error bars represent the dispersion (set 
to when only one template per spectral bin is available). 



10 <x depth of the 540s i' VW exposures is i' AB ~23.95, and 
we complement them by 420s z' exposures, with typical 10<x 
depths of z' AB ~22.8. Both sets of images are coadditions of 3 
subexposures separated by at least one night. We therefore have 
enough information to reject the vast majority of bad pixels, 
cosmic ray hits, and moving solar system objects. The time span 
of the 3 z' exposures on the other hand is usually too short to 
reliably recognize supernovae, which vary on time scales of a 
few weeks. We therefore again need to reject these contaminants 
at a later stage. 

Red-sequence Cluster Survey 2 (RCS-2) The RCS-2, 
designed to identify distant galaxy clusters through their 
galaxies on the red sequence, (lYee et al.l 12007b is an ongoing 
g'r'z' survey of 800 deg 2 at high galactic latitude to lower 
the absorption (see table Q] We have to date used ~600 deg 2 
kindly made available to us by the RCS-2 consortium, , and 
we complement their 360s z' band images by 500s or 680 
s i' exposures, depending on the seeing. The resulting 10<x 
depths of (z' AB ~22.5 and i' AB ~24.0 are similar to those of the 
CFHTLS-VW. Both the RCS-2 images and our complementary 
i' data are single exposures, which maximize the depth achieved 
for a given observing time. We use the RCS-2 g' and r' 
images, which are contemporaneous with the z' ones, to iden- 
tify and reject both supernovae and moving solar system objects. 

All images ar e pre-processed by the CFHT staff using the 
ELIXIR package dMagnier & Cuillandrdl2004T) . The CFHTLS 
Deep and Wide images are al igned and coadded by the Terapix 
data center dBertin et al.ll2002l) . For the CFHTLS Very Wide and 
RCS-2 datasets, we carry out our own processing to check and 
refine the astrometry and (for fields which overlap the SDSS) 
photometry. For the CFHTLS Very Wide, each pointing has 3 
subexposures per filter which are combined whilst rejecting bad 
pixels and cosmic ray impacts. The CFHTLS Very Wide and 
RCS-2 images in different filters are aligned (with distortion cor- 
rection) and trimmed to their common area. 

To date most of our volume coverage originates in the 
CFHTLS Very Wide and RCS-2, due in part to the late start 



of the z' part of CFHTLS-Wide. We have currently analysed 
350 deg 2 of the 800 deg 2 expected for these two surveys. Their 
relative shallowness has the advantage of producing targets for 
which spectroscopy can be obtained relatively easily on 8m-class 
telescopes. 

TableQ]summarizes the properties of the four surveys, listing 
the limiting magnitude, the maximum distances at which mid- 
L and late-T dwarfs can be detected to these magnitudes, and 
the current and final areas covered by the survey. Our full sur- 
vey probes several times the SDSS volume for T dwarfs and we 
expect to detect ~100 new T dwarfs (compared with the -150 
currently known) 

2.4. Near Infrared Imaging 

As explained above, we use /-band photometry to distinguish 
between brown dwarfs and z> 5.8 quasars. For brown dwarfs 
(and very low mass stars) the z' - J colour also provide a good 
spectral type diagnostic, for which we obtain better S/N ratio 
than i' - z' ■ That is in particular very helpful in eliminating mid- 
M dwarfs scattered into our i' - z' selection box by several <x 
noise excursions. Given the relative numbers of mid-M stars and 
brown dwarfs in a magnitude-limited sample, these noise excur- 
sions are sufficiently frequent to very significantly contaminate 
our i' - z' selection, but the better S/N ratio of the z' — J colour 
(and the low likelihood of large noise excursions at both i' - z' 
and z' - /) makes them obvious once we obtain /-band images. 

The /-band follow up has been carried out at several obser- 
vatories: La Silla (NTT, 3.6m), McDonald (2.7m), Kitt Peak 
(2.1m), La Palma (NOT, 2.5m), see|2] We adjust our integration 
times to achieve either a clear detection or a limiting magnitude 
that excludes any dwarf and demonstrates that the candidate is 
a high redshift quasar (usually around Jab = 22). For the few 
candidates which are not detected at i' and where we cannot 
exclude that the object was a supernovae in the z' image, we 
integrate deeper to detect the quasar at /. Any supernova has 
long faded, and cannot be detected at any reasonable depth. 

At the NTT, which accounts for most of our /-band follow 
up, we usually need exposures times of 5 to 10 minutes, obtained 
as ~ 40 seconds individual exposures which we jitter to measure 
and substract the sky background. The typical integration time 
on 2 meter-class telescope was 30 minutes. 

3. Candidate Selection 

Since brown dwarfs and high-redshift quasars are extremely red, 
and since our i' images are only moderately deeper than the z' 
ones, some of the most interesting targets are only detected at 
z' (i.e. are i'-dropouts). Any unrecognised artefact (cosmic ray 
impacts, unflagged hot pixel, optical ghost, etc..) in a z' image 
therefore translates into an invalid brown dwarf/quasar candi- 
date, since it associates a z' detection with (usually) an i' upper 
limit. Since brown dwarfs are very rare (~ 1 per 1 deg 2 MegaCam 
image, which contains from 50 000 to over 300 000 astronomi- 
cal sources), false detection rates of even 1 per 10 4 real sources 
would swamp true brown dwarfs in our candidate lists. We there- 
fore need to very effectively reject artefacts. Most of these are 
visually obvious, and we do examine every candidate before 
scheduling any follow-up observations, but the many hundred 
340 Megapixel images which we analyse contain too many arte- 
facts for this to be a practical first line of defense. 

Fortunately, high-redshift quasars (at the resolution of 
the MegaCam images) and brown dwarfs are point-like. We 
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Table 2. Technical characteristics of the telescopes used for the /-band follow up. 



Telescope Name 


Di timet er(m) 


In strument 


Exp time for 5 o~ de- 
tection of J(AB) = 
21 object 


Pixpl sralp 


Field of view 


NTT 


3.6 


SOFI 


45 s 


0.288" 


4.9' x 4.9' 


Kitt Peak 2.1m 


2.1 


SQIID 


~320s 


0.69" 


5.1' x5.3' 


Harlan-Smith 


2.7 


ROKCAM 


~3000s 


0.35" 


1.5' x 1.5' 


Telescope 












NOT 


2.5 


NOTCam 


~100s 


0.24" 


4.0' x4.0' 



therefore only need to distinguish point-like sources from both 
artefacts and extended objects, and don't have to tackle the much 
more difficult task of separating general astronomical sources 
from artefacts. Point Spread Function (hereafter PSF) fitting 
provides an excellent stellarity diagnostic, as well as optimal 
photometry and astrometry for point sources. It therefore forms 
the basis of our selection procedure. 



3. 1 . Image analysis 

We use the well known SExtractor (Bert in & Arnoutsl fl996) 
photometry package, to which two of us (Bertin & Delorme, 
in preparation) recently added a PSF-fitting module. In keeping 
with the general SExtractor philosophy, this module implements 
a dual-image mode, where source positions in a 'detection im- 
age' precisely determine where photometry will be extracted in 
a 'photometry image'. This dual image mode is particularly well 
matched to the extreme colours of our targets: given the relative 
depths of the i' and z' images, any object of interest is very ro- 
bustly detected at z' but faintly, if at all, at if. We therefore use 
the z' image as the detection image for both i' and z', naturally 
producing matched catalogues of i' and z' photometry for every 
object that is well detected at z', independently of its i' signifi- 
cance. This eliminates the delicate task of handling unmatched 
sources in independent catalogues: those might be weakly de- 
tected in the i' image, though with too low a significance for 
inclusion in any modestly reliable single image i' catalogue, and 
they therefore cannot validly be handled as pure upper limits. 

SExtractor implements simultaneous fitting of multiple PSFs 
to blended objects, providing accurate parameters for close bi- 
naries and usable measurements for point-like sources blended 
with galaxies. The latter is particulary important for the quasar 
search, since it can recover some lensed quasars which would 
otherwise be lost to confusion with their lensing galaxy. In ad- 
dition to more accurate parameters for the affected objects, this 
better handling of blends produces more complete catalogues in 
crowded fields. Introducing multiple PSFs recovers ~ 3% addi- 
tional sources in the relatively shallow CFHTLS-VW and RCS-2 
images, and ~ 10% in the deeper CFHTLS-Deep images. PSF- 
fitting also improves the photometric precision by ~10% over 
optimum aperture photometry (Bertin & Delorme, in prepara- 
tion), and it therefore allows us to use slightly deeper catalogues, 
for another 15% gain in sample size. Since low significance i' 
detections provide colours with complex error distributions, we 
replace them by the 5 cr detection limit on their image and com- 
pute a lower limit on if -z' ■ We note that the resampling involved 
in the coaddition of images built from multiple exposures, and 
in the filter to filter alignment, generates noise correlations on 
scales of 1-2 pixels. Thanks to the generous sampling of our 
MegaCam images (0.186"/pixel and seeing mostly above 0.6") 



and the use of a Lanczos 3 interpolation fun ction, source profiles 
are negligibly affected (Ber tin et al.ll2.002D . but resampling has 
a measurable low-pass filtering effect on photon noise. As we 
decided for practical reasons to ignore noise covariances in our 
fitting, the net effect on photometry is that errors estimates must 
be multiplied by a factor ~1.4. That factor is well determined 
for the CFHTLS-Deep and CFHTLS-Wide images, which are 
built from a large number of individual exposures, but for the 
CFHTLS-VW images it significantly varies from field to field 
according to the sub-pixel relative positions of the 3 coadded 
exposures. 

3.2. Filtering and target selection 

We start by requiring a > 10cr detection in the z' filter and a 
i' - z' > 1.7. These criteria without any additional filtering typ- 
ically yield over 10000 candidates per RCS-2 1 square degree 
field. The single exposures per filter used for the RCS-2 sur- 
vey (and for our follow-up of its fields) are most affected by 
cosmic ray hits and bad pixels, and the stacked images used in 
the other components contain fewer such artefacts. Flagging of 
known bad-pixel positions and simple morphological rejection 
of cosmic ray hits lowers this number to under 1000, but not to 
a point where visual examination would be practical. 

We then assess the stellarity of each candidate from its 
SExtractor output parameters to further decrease the number 
of false detections. After experimenting with several parameter 
combination, we have converged to the quality of the fit between 
the image and the PSF model, as summarised by the x 1 of the 
residuals, as our main diagnostic. We found that Sextractor's de- 
fault stellarity index, based on a specifically trained neural net- 
work, works well at high signal to noise ratios, but that it be- 
comes ineffective for the faint objects which dominate our cat- 
alogues. We similarly found that tests based on comparisons of 
fluxes through different apertures, or on peak surface brightness 
versus flux, are always less distinctive than^ 2 filtering. Since we 
currently prefer to visually inspect all final candidates, we very 
conservatively set our filtering threshold to a level where ~10% 
of the candidates are visually acceptable. These filtering criteria 
typically yield under 50 candidates per square degree. A lower 
threshold would not very significantly decrease the inspection 
workload, and might conceivably eliminate a few valid candi- 
dates. We will probably revisit this tuning as we gain experience 
with, and confidence in, our selection process. 

With our current settings, x 1 filtering reduces the number of 
artefacts by a factor of 10- 15 (figure. |5j- Visual inspection of the 
more than 2000 sources with i' - z' > 1 .7 in a 4 square degrees 
test region showed that^ 2 filtering rejected none of the 14 valid 
point-like sources, and our resolution of the few initial discrep- 
ancies was always in favour of the^ 2 filtering diagnostic. Further 
tests also verified that our current threshold is comfortably above 



P. Delorme et al.: Finding ultracool brown dwarfs with MegaCam on CFHT: method and first results 



7 



the highest^ 2 measured for valid candidates, and therefore very 
conservative. 

After this pruning of the initial catalogue to just bona fide 
point sources, we select candidates with an i' - z' criterion. 
Our synthetic photometry (figure [3) demonstrates that the L 
dwarf domain begins at i' — zf > 1.45. M dwarfs however, with 
i' - z' just below 1.45, massively outnumber brown dwarfs 
in a magnitude-limited sample. Poissonian photometric errors 
consequently scatter a significant number of M dwarfs into 
this L dwarf box. We therefore set our colour filtering to 
V - zf > 1 .7, or nominally to later than L4. With M8 and M9 
stars having i' - z' ~ 1.35, and assuming gaussian noise at our 
z' S/N limit, this colour threshold eliminates ~ 95% of these 
very late M dwarfs. Mid-M dwarfs have i' - z' ~ 1.1 and we 
eliminate ~ 99.9% of them. Since brown dwarfs are intrinsically 
much rarer, our candidate list nonetheless has some significant 
contamination by late-M dwarfs, but at a level which no longer 
overwhelms our follow-up capacity. 



4. Results 

4.1. Candidates 

We have so far analysed images for 357 deg 2 , in which we have 
identified 770 brown dwarf and quasar candidates. We have so 
far extracted /-band photometry for 2 1 5 of those, prioritising the 
reddest candidates (figure[6]l. 




Chi2 



Fig. 5. Histogram of the^ 2 of the psf fitting residuals. The thick 
red histogram shows all i' — z' > 1 .7 candidates after basic fil- 
tering, the medium-thick blue histogram those that additionally 
survived^ 2 filtering, and the thin green one the final candidates 
retained after visual inspection and supernovae rejection. These 
histograms show candidates from a 70 deg 2 test area. 



Our final filtering step is to eliminate as many as we can of 
the candidates which actually are supernovae that were bright at 
the time of the z' image. All RCS-2 fields have contemporane- 
ous g' , r' and z' images, and for supernovae the g' and r' im- 
ages are very significantly deeper than the z' one. We therefore 
very reliably reject their supernovae by inspecting these g' and r' 
images. The CFHTLS-Deep images are stacks of exposures ob- 
tained over several years, and any supernova is eliminated by the 
sigma-clipping applied during stacking. The exposures which 
contribute to our CFHTLS-Wide and CFHTLS-VW z' images, 
on the other hand, were usually obtained over shorter time spans 
than the ~6 weeks dPain & SNLS Collaborationl2003l) timescale 
of supernovae photometric evolution. We therefore mostly can- 
not recognize their supernovae based on their photometric vari- 
ation between the individual exposures, and the exposures in 
other filters are usually not sufficiently contemporaneous to re- 
ject them based on a blue instaneous colour. We must therefore 
handle some supernovae contamination at a later stage. 



150 



100 




Fig. 6. Histograms of the number of if - z' candidates as a func- 
tion of i' - z' colour (light green), and of the number of these 
candidates for which /-band photometry is currently available 
(dark red). 



This first set of i'z'J photometry allows a good assessment 
of the actual nature of these 215 candidates. As summarised in 
Table [3] they include 23 likely T dwarfs, 57 L dwarfs and very 
late M dwarfs candidates (M8 and M9 dwarfs have very similar 
z' - J colours to L dwarfs), and at least 4 high redshift quasars 
(published in IWillott et al.ll2007l) . 22 targets remain undetected 
in deep /-band images and are most likely artefacts which our 
filtering did not catch. 109 objects have z' - J < 1.6 and most 
of those are likely M dwarf contaminants. Some however have 
/-band upper limits which are insufficiently deep to ascertain 
whether they are artefacts, quasars or Mid-M dwarfs. Those will 
need additional follow-up to clarify their status. 

Since we prioritised the analysis of the reddest candidates, 
this first i'z'J sample is strongly biased towards T dwarfs. 
Appproximately correcting for this bias, we estimate that 
~40% of our i' - z' candidates are actual cool dwarfs, of 
which ~ 15% are T dwarfs. The T dwarfs include several with 
extreme colours, which ongoing spectroscopic observations will 
characterize further. One, CFBDS J003402-005206, actually 
is an independent discov ery of ULAS J003402.77-005206.7 
which IWarren et al.1 d2007l) recen tly identified with the UKIDSS 
survey dBurningham et al] 120071) as a brown dwarf later than 
T8 (~ 65QK). Even more recently, we obtained spectro- 
scopic observations of an even cooler brown dwarf (~ 625K, 
lDelorme"e t al. 2008) Figure Q shows several other candidates 
that are at least as promising, and those are currently queued for 
near-IR spectroscopy. 
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Table 3. Preliminary classification of the candidates with /-band photometry 





T dwarfs 


L dwarfs 


quasars 


M dwarfs 


artefacts 


Number of objects 


23 


57 


4 


109 


22 


Percentage 


10.7% 


26.5% 


1.9% 


50.7% 


10.2% 



This classification is based on their position in the i' - z' versus z' — J colour-colour diagram. 



Y dwarfs ?- 




IV z" colour (AB mdy) 

Fig. 7. i' - z' vs z' - J colour-colour diagram of our cool dwarf 
candidates (black). The error bars are 1 cr when we have i' 
and z' measurements. An arrow is used when i' - z' is only 
a lower limit. For clarity, the z' - J errors are not shown; 
they are usually under 0.15 mag. The spectra of the high- 
lighted T3, T3.5 and T5 objects are presented in this article. 
The big blue squares represent the coolest brown dwarfs known, 
ULAS J0034 02.77-0052067/CFB DS J003402-005206, recently 
published by Warr en et al.l ((2007) and independently identified 
in our survey and CF BDS J005910.83-01 1401.3, presented in 
iDelorme et al. (120081). The T7 object highlighted here was pub- 
lished by IChiu et alX(l2006l) while queued for our CFBDS spec- 
troscopy. The spectral type domains are based on our synthetic 
photometry. 



tion correction. The spectra were extracted and calibrated using 
our own IDL procedures. 



The reduction proceeds as follows. The sequence of spectral 
images are flat-fielded using an internal flat taken immediately 
after the science frames. The five useful cross-dispersed orders 
are then extracted in five individual images that are corrected for 
distortion in the spectral dimension. For most objects, the trace 
is too faint over many wavelengths intervals to determine trace 
position, so its curvature is derived from the reference star spec- 
trum. These individual order frames are then pair-subtracted, ef- 
fectively removing most of the sky, dark current and hot pixels 
contributions. Each frame is then collapsed along the spectral di- 
mension to determine the positive and negative traces positions. 
We then extract the spectra using positive and negative extraction 
boxes that have identical but opposite integrals; this minimizes 
the contribution from residual sky line that would remain from 
the pair subtraction. The same operation is performed on the AO 
telluric calibration star. Spectra derived from individual image 
pairs are then median-combined into final target and calibration 
star spectra. A telluric absorption spectrum is then derived using 
the calibration star spectra. A black body spectrum with a tem- 
perature of 10 000 K is assumed for the AO stars and hydrogen- 
lines are interpolated over. The target spectrum is then divided 
by the derived telluric transmission spectrum. A first order wave- 
length calibration is obtained from an argon-lamp spectrum, and 
fine-tuned by registering the bright OH-lines obtained from a 
sum of the pair of images of interest. Table |4]summmarizes the 
properties of the 3 objects. 



4.2. Spectroscopy 

We present here near-IR spectroscopic observations of our three 
first T dwarf candidates, as an illustration of the content of our 
candidates catalogu e. Two, which we originally published in 
IWillott et all J2005), originate in the CFHTLS-Deep survey and 
are thus fainter than most of our candidates. The last comes 
from the Very Wide component of the CFHTLS and is more 
representative. Cross dispers ed spectra were obtained during 
semester 2006 A with GNIRS (Eli aset al.l2006l) at Gemini South 
in Service mode. The slit width of 0.68 arcsec coupled with 
the short camera and the 31.7 1/mm grating yielded a resolv- 
ing power of 900, and the spectra have full wavelength coverage 
between 0.9 and 2.4 microns. A-B (not ABBA) sequences were 
used, with individual 5 -minutes exposures for the brighter target 
and 10-minutes for the fainter CFHTLS-Deep targets. The total 
on-source integration time is 30 minutes for CFBDS 193430 - 
214221 (Jab = 17.9), 180 minutes for CFBDS 1001 13 + 022622 
(Jab = 19.7) and 200 minutes for CFBDS0959U + 023655 
(Jab - 20.3). The OH sky lines were used for wavelength cal- 
ibration, and an A-type star was observed immediately before 
each sequence for relative flux calibration and telluric absorp- 



The spectra (figure. [8) confirm that all three candidates with 
spectroscopic observati ons are T-dwarfs. We determined their 
spectral types using the Burgasser et alj d2006f) spectral indices. 
Table[5]lists thoses indices and the corresponding spectral type. 
We retain as our prefered determination the "weighted" spec- 
tral types, rounded to the closest half-integer. These "weighted" 
spectral types take into account the better sensitivity of those 
indices that vary most for a given subtype range. The reddest 
target, CFBDS 1001 13+022622 (? - z' > 4.2), turn out to be 
a T5±l dwarf, and its 45 to 110 parsecs photometric distance 
makes it one of the farthest mid/late-T dwarf currently known. 
The large uncertainty on its distance is dominated by the spec- 
tral typ e uncertainty and the >1 dimming between T4 and T6 
dwarfs. dVrba et al] 12004). with photometric uncertainties con- 
tributing less than 5%. However, the absolute magnitudes of 
mid-T dwarfs are themse lves uncertain b y as much as 1 mag- 
nitude (see for instance iLiu et alj 120061) . The faintest of the 
three dwarfs, CFBDS0959 14+023655 (f - z' > 3.6) turns out 
to have an earlier spectral type, T3, but lies even farther, be- 
tween 120 and 130 parsecs. The indices of the brighest one, 
CFBDS 193430-214221 (i'-z' > 3.4), indicate a T3 .5+0.5 spec- 
tral type. The spectral type uncertainties are derived from the 
scatter between the estimates from the various spectral indices. 
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5. Conclusion 

Our survey has to date found 23 T dwarf candidates, 57 L or 
very-late M dwarf candidates, and 4 high redshift quasars, out of 
215 candidates with i', z' and J magnitudes. These were drawn 
from a larger sample of 770 candidates with i' and z' magni- 
tudes, found in 357 deg 2 of i' and z' MegaCam images. Taking 
into account our prioritising of the reddest candidates for /-band 
observations, we expect that complete follow-up of these 770 
candidates will yield ~45 T dwarfs and 200 L dwarfs. Scaling 
to our final ~ 800deg 2 of shallow surveys, RCS-2+Very Wide, 
then predicts ~ 100 T dwarfs and over 450 L or very late-M 
dwarfs, approximately doubling the number of known brown 
dwarfs. Our analysis of the CFHTLS-Deep and CFHTLS-Wide 
surveys has, and will, yield additional candidates at large dis- 
tances, which will constrain the galactic scale height of brown 
dwarfs. We plan to obtain spectra for the most exciting of these 
many brown dwarfs, and expect that the large discovery volume 
will produce eve n cooler objects than our recent T9/Y0 discov- 
ery, described in lDelorme et alj d2008l) . 

Appendix A: Finding Charts 
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Name 




RA 


DEC 


Jab 


z' AB - Jab 


'as z ab 


Spectral 


Type 


Distance Range(pcs) 


CFBDS 100 113+022622 


10 


01:13.1 


+02:26:22.3 


19.7 


3.1 


3.75 


T5± 


1 


45-110 


CFBDS 193430-214221 


19 


34:30.4 


-21:42:21.0 


17.9 


2.8 


> 3.4 


T3.5d 


:.5 


38-48 


CFBDS0959 14+023655 


09 


59:14.8 


+02:36:55.2 


20.3 


3.1 


3.7 


T3.0d 


:.5 


120-130 



Table 5. Spectral indices and the resulting spectral types 



Spectral Index 


CFBDS 1001 13+022622 


CFBDS 193430-214221 


CFBDS0959 14+023655 


CH 4 -J 


0.43 - T5.2 


0.59 - T2.9 


0.57 - T3.2 


H 2 0-/f 


0.46 - T3.0 


0.48 - T2.7 


0.45-T3.1 


CH 4 -H 


0.40 - T5.7 


0.63 - T3.9 


0.79 - T3.0 


CH 4 - K 


0.27 - T4.7 


0.34-T4.1 


0.54 - T3.0 


Average Spectral Type 


T4.6±1.2 


T3 .4+0.7 


T3. 1+0.1 


Weighted Spectral Type 


T5.0 


T3.6 


T3.1 


Final Spectral Type 


T5.0±l 


T3.5+0.5 


T3.0+0.5 



Spectral types are derived from Burg asser et al.l (|2006) spectral indices. 
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Fig. 8. GNIRS near-IR spectra of the three T Dwarfs observed on Gemini-South. The two fainter targets originat e from the CFHTLS- 
Deep, while the brighter ones comes from the CFHTLS-VW. Templates of spectral types T2 and T4.5 from lMcLean et aD d2003h 
are shown for comparison. All spectra are normalised at 1.26 microns and displayed with integer offsets for clarity. 



